Electronic structure and Fermi surface character of LaONiP from first principles 
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Based on First-principles calculation, we have investigated electronic structure of a ZrCuSiAs structured su- 

" QQ ' perconductor LaONiP. The density of states, band structures and Fermi surfaces have been given in detail. Our 

^^ , results indicate that the bonding of the La-0 and Ni-P is strongly covalent whereas binding property between 

^^ ■ the LaO and NiP blocks is mostly ionic. It's also found that four bands are across the Fermi level and the corre- 

^vq ' sponding Fermi surfaces all have a two-dimensional character. In addition, we also give the band decomposed 

charge density, which suggests that orbital components of Fermi surfaces are more complicated than cuprate 

superconductors. 
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PACS numbers: 71.18.+y, 71.15.Mb,74.25.Jb 



I. INTRODUCTION 



i^ . Quaternary rare earth transition metal phosphide oxides with the tetragonal ZrCuSiAs structure^ are known for more than 

Q ■ twenty years^. Numerous quaternary phosphide oxides have been reported, including the actinoidand copper-containing com- 
I \ pounds UOCuP,^ ThOCuP'^as well as the series of lanthanoid (Ln) compounds LnOFeP^, LnORuP^, LnOCoP^, LnOOsP, 
P^- LnOZnp5'6 and LnOMnP^ . 

^ \ However, most of interests are limited to the synthesization and structural analysis, only few physical properties have been 

^ . investigated for this compound series. Because of open- shell structure and electronic correlation of 3d electrons, the ZrCuSiAs 
type transition metal phosphide oxides are expected to display interesting magnetic and electronic phenomena. With the same 
crystal structure, the different groundstates including antiferromagnetic insulator^, semiconductor and metal^have been found 
in the ZrCuSiAs type compounds. Moreover, LaOFeP^'^^ has been found to be a type II superconductor with Tc of 3.2~6.5 K , 
F^ \ which leads to renewed interest in studying physical properties in this class of materials. Later, Tc around 26 K has been found 
^ ■ in F-doped^^and Sr-doped^^ LaOFeAs. Very recently, up to 7^^=55 K was reached in LnOi-j^FxPeAs through replacement of 
O I La by Ce^^, Pr^"^, Nd^^, Sm^^'^^,and Gd^^. Ferromagnetic Kondo lattice systems found in Ce-based phosphide oxides CeOTP 
^ O^ ■ (T=Ru,Os)i2fe20 also attract much attentions. 

Recently, a new compound LaONiP with ZrCuSiAs structure has been reported by Watanabe et al?^, and it has also been 
^s^ . synthesized by Tegel et al^^ with different experimental method. It has been identified as a superconductor with a critical 
^ ■ temperature Tc of 3 and 4.3 K, respectively. Superconductivity was also found in LaONiAs and LaOi-j^F^^NiAs^. LaONiP 
OO , is the first nickel phosphide oxide and the second superconducting compound in LaOMP(M=transition metal). It seems that 
C/N ■ the superconductivity may be a common phenomena in this family. Although this Tc is lower than that of the copper oxides, 
because of possible different underlying mechanism, the discovery of superconductivity in the new material systems can provide 
valuable knowledge for understanding of superconduction and for finding another superconductor. Thus, it is crucial to clarify 
"^" ' the underlying mechanism of novel superconductivity in this series of compounds. Whereas the understanding of electronic 
structure especially the Fermi surface character is a basic step to obtain insight into the superconductivity mechanism. 

In present calculation, we perform the systematic investigation for LaONiP based on density functional theory. The de- 
tailed structural parameters, electronic structures and Fermi surfaces were given. The bonding properties are also discussed in 
combination with the concept of two-dimensional (2D) building blocks and a Bader analysis of the charge density. The band 
decomposed charge density is also provided to analyze the orbital character of each Fermi surface. 



II. COMPUTATIONAL DETAILS 

The present calculations have been performed using the Vienna Ab initio Simulation Package (VASP) code^"^'^^ within pro- 
jector augmented-wave (PAW) method,^^'^^ general gradient approximations (GGA)^^were used in the present calculations. The 
\.di(5s^5p^5d^6s^), Ni(3J^4^^), P(3^^3/7^), Oils^lp"^) are treated as valence states. To ensure enough convergence, the energy 
cutoff was chosen to be 600 eV, and the Brillouin zone was sampled with a mesh of 16x16x8^ points generated by the scheme 
of Monkhorst-Pack^^. A first order Methfessel-Paxton method with cr = 02eV was used for relaxation. The crystal cell and 
internal parameters were optimized using the conjugate gradient method until the total forces on each ion less than 0.01 eV/A. 
Then density of states (DOS) calculations were performed using the tetrahedron method with the Blochl corrections^^. The 
Fermi surfaces and 3D charge density iso-surfaces have been drawn by Xcrysden^. 
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FIG. 1 : (Color online)Crystal structure of LaONiP and near-neighbour coordinations of La and Ni atom. The element symbols have been 
labeled in the figure. 
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III. RESULTS AND DISCUSSION 



A. Crystal Structure 



LaONiP crystallizes in the tetragonal ZrCuSiAs structure (space group P4/nmm) with the Ni atom at 2a (0.75, 0.25, 0), O at 
2b (0.75, 0.25, 0.5), La at 2c (0.25, 0.25, zi) and P at 2c(0.25, 0.25, zi)- It has a layered structure composed of an alternating 
stack of LaO and NiP block. The crystal structure is therefore rather simple with eight atoms (two formula units) in the unit 
cell as shown in Fig.[T] Lanthanum is eightfold coordinated by four oxygen and four phosphorus atoms. The nickel atoms are 
tetrahedrally coordinated by four phosphorus atoms, forming a distorted tetrahedra with two different P-Ni-P anglesi^i^^^ And 
they also have four neighboring nickel atoms within the same layer. 

Table [J gives the optimal structural parameter of LaONiP together with available experimental results. The experimental 
results are very consistent with each other. The present calculation is also in agreement with experiment, which indicates 
that the calculation is accurate. The slight difference between experiment and calculation may be induced by the fact that the 
LDA/GGA calculation cannot describe the electronic correlation between 3d electrons of Ni atom accurately. It should be 
noticed that the electronic structure ^^^and magnetic properties'^ in case of LaOFeAs are very sensitive to small change of 
internal coordinates, as shown in Ref . I32ll33ll34l . It seems that in order to give a more reasonable description, further calculation 
including the electronic correlation should be performed. 



B. Density of States and Binding Properties 

Fig. [2] shows the calculated total and partial DOS of LaONiP. We can see the states below -13 eV is contributed by La and O 
atom. And the La-p and 0-s shows strong hybridization with each other under the range -20 to -13 eV, indicating strong La-0 
bonding. The energy range from -12 to -10 eV is predominated by the P and Ni atoms with only a weak mixture of La atom. 
Whereas all atoms are contributing to the DOS at the energy range -6 to -3 eV. It should be noticed that from -2 to 2 eV, the 3d 
states of Ni atoms dominate the DOS, together with a slight contribution from P-3p states, while the La and O almost may be 
neglected. These results indicate that there are a significant overlap between the orbit in La-0 and Ni-P and the bonding differs 
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considerably from ionic ones. 

Recently, the electronic structure of LaOFePn(Pn=P,As)^2i^^i2i2^have been investigated extensively. These calculations 
predict that the Fermi energy lies just above a peak in the DOS. Thus, LaOFePn have a very steep and negative slope of the 
density of states (DOS) at the Fermi level, which drives the system close to a magnetic instability. Compared with iron-based 
compound, we find that the shape of DOS in case of LaONiP is similar. However, because of Ni^^ (3d^) contributes two more 
electrons than Fe^^ (3d^)in LaOFePn, Fermi level is lifted up in case of LaONiP and the densities of states don't decrease 
monotonously with energy, which may put LaONiP away from magnetic instability. 

As shown in above DOS analysis, we can find that the states of Ni and P atom predominate DOS near Fermi level, which 
suggests that the superconductivity origins from the contribution of the Ni 3d and P 3p as shown in Fig. |2(a)[ The layered 
copper-based superconductors have been extensively studied for several decades. In cuprate superconductors, Cu^^ occupies a 
planar 4-fold square site and the charge carriers at the Fermi level are driven by the d^^-yi orbit. While Ni^^ in LaONiP occupies 
a tetrahedral site coordinated with four P^" ions. Such a marked diff'erence in the coordination structure between the LaONiP and 
cuprate superconductor is expected to lead to the diff'erent mechanism of superconductivity. The Ni ion in LaONiP is formally 
in a 3 J^ configuration. In a tetrahedral crystal field. The Ni d bands will split into a lower lying e and an upper lying t2 orbit. 
However, because of the distorted tetrahedral as shown in Fig. |2(b)[ a clear separation in energy of the d orbit in a lower e Od^i, 
3dx2-y2) and higher t2(3dxy, 3dyz, 3dxz) set is not seen in the present calculation. Whereas 3dxy and 3d^2 contributes more to DOS 
near the Fermi surface, which may play a more important role in superconductivity. In addition, from Fig. |2(b)^ nd Fig. |2(c)[ all 
five 3d orbits of Ni atom and three p orbits of P atom has comparable contribution to DOS at the Fermi level. 

We also analysis the DOS using the concept of two-dimensional (2D) building blocks^^'^^'^^"^^. It will be helpful to understand 
how charge transfer and hybridization eff'ects will aff'ect their electronic and binding properties when building the solid from 
these separated blocks. From the Fig. [3l we can see that the LaO states are pushed up in energy by approximately 2 eV , and 
becoming insulating in the solid but the isolated LaO layer is metallic, which indicate there are electron transfers from the LaO 
to NiP when building the solid. From the figure, we also conclude that the LaO-FeP interaction is strongly ionic character with 
weak hybridization under energy range of -5 and -2 eV. 

As shown above, the qualitative electronic structure and binding properties is very similar to the results of LaOFeP give in 
Ref. [35I. In order to understand the binding properties more quantificationally, we also performed a Bader analysis"^^"^^ of the 
charge density. The core charge missed in general pseudopotentials method has also been added to perform this analysis. Table|IIl 
gives the charges of each atom specie using bader analysis together with the pure ionic picture. We can notice that the La-0 and 
Ni-P bonds diff'er considerably from ionic character, whereas the bonding between the LaO and FeP blocks is mostly ionic. It 
is consistent with the previous DOS analysis. Compared with the result of LaOFeP^^, we also find that the deviation of ionic 
character between Ni and P is larger than the case of LaOFeP, which suggests the co valency between 3d metal and P atoms is 
stronger in LaONiP. The weak p-d hybridization between Fe and As is also found in case of LaOFeAs^SiSi^. 

C. Band Structure and Fermi Surface 

Now, Let's turn our attention to band structure. As shown in Fig.lH the characteristic feature of band structures is the strongly 
pronounced two dimensionality along F-Z and A-M. We observe four bands across the Fermi energy, which have been indicated 
in the figure with diff'erent colors. The corresponding Fermi surface in the first Brillouin zone are displayed in Fig. \5\ Due 
to the two-dimensional electronic structure, all Fermi surfaces are cylindrical-like sheets parallel to the k^ direction. The first 
sheet is hole cylinder centered along R-X whereas the other three are electron cylinders centered along the A-M high symmetry 
line. Such a series of Fermi surfaces are diff'erent from the results of LaOFeP^^ and LaOFeAs^^'^^'^^, where five bands across 
the Fermi level with four of them being cylindrical-like sheets(two electron cylinders along A-M and the other two hole-like 
cylinders along F-Z direction) and the other one being a three-dimensional hole pocket along F-Z. The absence of distorted 
sphere Fermi surface sheet suggests the two dimensionality of band structure in LaONiP is more pronounced than the case of 
LaOFePn. This can be explained by the fact that Ni^^ (3d^) contribute two more electrons than Fe^^ in LaOFePn , nickel based 
system has relatively higher carrier density and the hole bands tend to be filled. Thus, three-dimensional hole pocket disappeared 
and electron bands dominate the conductivity. This is also been evidenced by recent experiment^^, in which they found that the 
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FIG. 2: (Color online) Total density of states and Partial density of states for each atomic species(a) and the orbit resolved density of states for 
Ni 3d states(b) and P 3p states(c) of LaONiP. The Fermi level is at zero energy. 
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FIG. 3: (Color online) Comparison of the DOS of the LaO(left) and NiP (right) blocks calculated as isolated entities with the PDOS for LaO 
(center left) and NiP (center right) in the LaONiP solid. 




FIG. 4: (Color online) The calculated band structure of LaONiP. The Fermi level is at zero energy and marked by a horizontal dashed line. 
Four bands across the Fermi level are indicated in figure and marked with diff'erent colors. 



charge carriers of LaONiAs are dominantly electron- type. 

For further insight into the mechanism, we give the band decomposed charge density in Fig.O Although the Fermi surface 
of all bands have similar shape, the iso-surfaces of band decomposed charge density are different from each other as shown in 
Fig. [6] and Fig. [71 which also suggests that the orbital component is different. From the shape of iso-surface, we can speculate 
that the band a is made up of Ki-?>dxy and weakly mixture with p - cr orbit, which may be built with px and py. In case of the 
band J, the Ni orbit is mainly contributed by 3(i^2_^2 and the states of P atom exhibit a p^ orbital character. It seems that the 
orbital component of band b and c is more complicated than a and d. Combined with different orbital picture, we can deduce the 
3d orbital character of c is mostly contribute by 3d^2 and maybe some component of 3dx2-y2 . Whereas the band b may be mixed 
with 3dxz and 3dyz, the states of P are similar to a, which may be composed of px and py. Our results shows that not only the 
3dx2-y2 driving superconductor in layered copper-based oxides, but also all five 3d orbits are contribute to charge carriers at the 
Fermi level. So, it seems that the underlying mechanism is very complicated in ZrCuSiAs structured superconductor. 
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FIG. 5: (Color online) The Fermi surface sheet of LaONiP shown in the first Brillouin zone centered around the F point. 





FIG. 6: (Color online) The band decomposed charge density of four bands across the Fermi level. Iso surfaces correspond to 0.035 e/A^. In 
order to visualize the orbital character of P atoms, we choose 0.025 e/A^ for band a. 












FIG. 7: (Color 
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online) Top view of 
035 ^/A^ In order to 
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the iso-surface of band decomposed charge density of four bands across the Fermi level.Iso surfaces 
visualize the orbital character of P atoms, we choose 0.025 e/A^ for band a. 



IV. CONCLUSION 

Based on First-principles calculation, we have investigated the electronic structure, binding properties and Fermi surface 
character of a new ZrCuSiAs structured superconductor LaONiP. Our results indicate that the density of states and bonding 
properties are similar to the results of LaOFeP. However, the Fermi surface character exhibits remarkable difference and two- 
dimensional character is more significant in case of LaONiP. We also find that the orbital components of Fermi surfaces are 
very complicated, which suggests that the superconductivity mechanism is remarkably different from the well-known cuprate 
superconductors. 
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